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Abstract. We address a question about how much resonance poles and residues extracted from data depend on a model used 
for the extraction, and on the precision of data. We focus on the Py \ TtN scattering and use the dynamical coupled-channel 
(DCC) model developed in Excited Baryon Analysis Center (EBAC) at JLab. We examine the model-dependence of the poles 
by varying parameters largely within the EBAC-DCC model. We find that two poles associated with the Roper resonance are 
fairly stable against the variation. We also study the stability of the Roper poles against different analytic structure of the P\ \ 
amplitude below nN threshold by using a bare nucleon model. We again find a good stability of the Roper poles. 
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INTRODUCTION 

Extraction of N* information, such as pole positions and vertex form factors, is an important task in hadron physics. 
This is because they are necessary information to address a question whether we can understand baryon resonances 
within QCD. In order to extract the N* information, first, one needs to construct a reaction model through a com- 
prehensive analysis of data. Then, pole positions and vertex form factors are extracted from the model with the use 
of the analytic continuation. Therefore, the N* information extracted in this manner is inevitably model-dependent. 
There are several different approaches to extract the N* information. Although almost all 4-stars nucleon resonances 
listed by Particle Data Group (PDG) are found in all approaches, existence of some N* states, in particular those in 
the higher mass region, is controversial. Thus, commonly asked questions are how much model-dependent the ex- 
tracted resonance parameters are, and how precise data have to be for a stable resonance extraction. These are the 
questions we address at Excited Baryon Analysis Center (EBAC) at JLab[l], within a dynamical coupled-channels 
model (EBAC-DCC) [2]. We focus on the nN P\ \ partial wave and the stability of its pole positions, particularly those 
corresponding to the Roper resonance. In the region near Roper A^(1440), two poles close to the nA threshold were 
found in our recent extraction [3] from the JLMS model [4] (JLMS is one of EBAC-DCC model), while only one pole 
in the similar energy region was reported in some other analyses. We examine the stability of this two-pole structure 
against the following variation, keeping a good reproduction of SAID single-energy (SAID-SES) solution [5] unless 
otherwise stated. 

• Large variation of the parameters of the meson-baryon and bare N* parameters of the EBAC-DCC model. 

• Inclusion of a bare nucleon state: The analytic structure of this model is rather different from the original EBAC- 
DCC model, in particular in the region near the nucleon pole [6], . 

• Fit to the solution based on the Carnegie-Mellon University-Berkeley model (CMB) [7] which has rather different 
behavior from SAID-SES for higher W. 

DYNAMICAL COUPLED-CHANNELS MODELS 

Here, we briefly describe dynamical coupled-channels models; the EBAC-DCC model and the bare nucleon model. 

The EBAC-DCC model contains nN, T]N and nnN channels and the nnN channel has nA, pN and oN components. 
These meson-baryon (MB) channels are connected with each other by meson-baryon interactions (vmb.m'b'X or 
excited to bare N* states by vertex interactions (Tmb^n*)- With these interactions, the partial-wave amplitude for 
the M{k) + B(-k) — > M'(k')+B'(-k') reaction can be written by the following form: 

TMB,M'B'{k,k' ,E) — t MBM i B i{k,k\E) + t^ BM , B ,(k,k',E), (1) 



where the first term is obtained by solving the following coupled-channels Lippmann-Schwinger equation: 

tMB,M'B'(k,k' ,E) = v M B,M'B'(k,k')+ / q 2 dqv M B,M"B"(k,q)G M ii B ii(q,E)t M it B i/ M i B i(q,k' ,E). (2) 

Here Cmb is the integration contour in the complex-^ plane used for the channel MB. The second term of Eq. (1) is 
associated with the bare N* states, and given by 

tMBtf'B'fcV > E ) = ^r MB ^ N *(k,E)[D(E)} u T N *^ M/B ,(k',E), (3) 

where the dressed vertex function F N *^ M r B i(k,E) is calculated by convoluting the bare vertex T N *^ M i B i{k) with the 
amplitudes t MB M i B i(k,k' ,E). The inverse of the propagator of dressed N* states in Eq. (3) is 

[D-\E)} iJ = (E-m Nt )S i j-L i j(E), (4) 
where is the bare mass of the i-th N* state, and the N* self-energy is defined by 

*^mb (q ! E)Gmb (q i E)Tm B ^n* (q,E). (5) 

mb Jc mb 

Defining E a (k) = [m 2 a + k 2 ] 1 ! 2 with m a being the mass of particle a, the meson-baryon propagators in the above 
equations are: GMB(k,E) = 1/[E— £m(£) — E B (k) + ie] for the stable nN and ijN channels, and GMB{k,E) = 
l/[E — EM(k) —E B (k)—T.MB(k,E)] for the unstable kA, pN, and aN channels. The self energy Lmb (k, E) is calculated 
from a vertex function defining the decay of the considered unstable particle in the presence of a spectator noxN with 
momentum k. For example, we have for the ftA state, 

y (h f\ mA ( ,2, M„ N (g) L/A-Wg)! 2 f „ 

L ^ t} E A (k)J Ci q dq [Ml N {q)+k 2 Yl 2 E-EAk)-[Ml N {q) + k 2 Yl 2 + ^ W 

where M nN (q) — E n (q) +E^(q) and /a->jtjv(?) defines the decay of the A — > nN in the rest frame of A, C3 is the 
corresponding integration contour in the complex-^ plane. 

To examine further the model dependence of resonance extractions, it is useful to also perform analysis using 
models with a bare nucleon, as developed in Ref. [8]. Within the formulation of EBAC-DCC model, such a model 
can be obtained by adding a bare nucleon (No) state with mass and — > MB vertices and removing the direct 
MB — > N — > M'B' in the meson-baryon interactions Vmb.m'b'- AH numerical procedures for this model are identical to 
that used for the EBAC-DCC model, except that the resulting amplitude must satisfy the nucleon pole condition: 

R , , , p \ [FnNN(kon)] 2 

'nN nNK^ ~ ' Komk ~~ * Kon^ ~^ m N) — ~~ a ~T7 t- (I) 

E-m%- L(wiv) 

with 

m N = m% + t(m N ) and F xNN (k oa ) = F^(k on ) . (8) 



Here we have used the on-shell momentum defined by E = y mjj + k 2 n + yjm^+k 2 ^ Also, is the self-energy 

for the nucleon. More details for the calculational procedure following Afnan and Pearce is found in Refs. [1, 8]. 



RESULTS 

Now we show our numerical results to examine the stability of the P\i poles. In the following subsections, we present 
results from various fits by varying the dynamical content of the EBAC-DCC model, and by using a model with a 
bare nucleon. We show in figures the quality of fits of these models, and in Table 1 the pole positions from the models 
as well as % 2 per data point (Xp d )- We find the poles with the method of analytic continuation discussed in detail in 
Refs. [3, 9]. In Table 1, we also present pole positions from JLMS[4] and SAID-EDS (energy-dependent) [5]. 



TABLE 1. The resonance pole positions Mr for P\ \ [listed as (ReM^, — IitiMr) in the 
unit of MeV] extracted from various parameter sets. The location of the pole is specified 
by, e.g., (sxNi s riNi s 7C7zNj s KAj s pNi s oN) = (upuupp), where p and u denote the physical 
and unphysical sheets for a given reaction channel, respectively. is x 2 P er data point. 



Model 


upuupp 


upuppp 


uuuupp 


uuuuup 




SAID-EDS 


(1359, 81) 


(1388, 83) 






2.94 


JLMS 


(1357, 76) 


(1364, 105) 




(1820, 248) 


3.55 


2N*-3p 


(1368, 82) 


(1375, 110) 




(1810, 82) 


3.28 


2N*-4p 


(1372, 80) 


(1385, 114) 


(1636, 67) 


(1960, 215) 


3.36 


liV ljV*-3p 


(1363, 81) 


(1377, 128) 




(1764, 137) 


2.51 




2A/*-3p and 2A/*-4p fits 

We varied both the parameters for the meson-baryon interactions (v MB M i B i) and parameters associated with bare N* 
states (mjL, rW*<_>Affl)- The obtained meson-baryon interactions are quite different from those of JLMS. We obtained 
several fits which are different in how the oscillatory behavior of SAID-SES amplitude for higher W is fitted. The 
results from the 2/V*-3p (dotted curves) and 2/V*-4p (dashed curves) fits are compared with the JLMS fit (solid curves) 
in Fig. 1. The resulting resonance poles are listed in the 3th and 4th rows of Table 1. Here we see again the first two 
poles near the 7rA threshold from both fits agree well with the JLMS fit. This seems to further support the conjecture 
that these two poles are mainly sensitive to the data below W ~ 1.5 GeV where the SAID-SES has rather small errors. 
However, the 2/V*-4p fit has one more pole at Mr = 1630 — /45 MeV. This is perhaps related to its oscillating structure 
near W ~ 1.6 GeV (dashed curves), as shown in the Figs. 1. On the other hand, this resonance pole could be fictitious 
since the fit 27V* -3p (dotted curve) with only three poles are equally acceptable within the fluctuating experimental 
errors. Our result suggests that it is important to have more accurate data in the high W region for a high precision 
resonance extraction. 



Here we show our results obtained with the bare nucleon model, and then address the question whether difference 
in the analytic structure of the nN amplitude below nN threshold strongly affects the resonance extractions. The bare 
nucleon model is fitted to SAID-SES, and at the same time, to the nucleon pole conditions Eqs. 8. Meanwhile, the 
original EBAC-DCC model has different singular structure below the nN threshold. The question is whether such 
differences can lead to very different resonance poles. Our fit of the bare nucleon model is shown in Fig. 2 and 




FIGURE 2. The real (left panel) and imaginary (right panel) parts of the P\ \ amplitudes. 



compared with SAID-SES and JLMS. We see that the two fits agree very well below W = 1.5 GeV, while their 
differences are significant in the high W region. The corresponding resonance poles are given in Table 1 . We also see 
here that the first two poles near the 7rA threshold are close to those of JLMS. Our results seem to indicate that these 
two poles are rather insensitive to the analytic structure of the amplitude in the region below nN threshold, and are 
mainly determined by the data in the region m# + m K < W < 1.6 GeV. 



CONCLUSION 

We have examined the stability of the two-pole structure of the Roper resonance. We showed that two resonance poles 
near the 7tA threshold are stable against large variations of parameters of meson-exchange mechanisms within EBAC- 
DCC model [2]. This two-pole structure is also obtained in an analysis based on a model with the bare nucleon state. 
Our results indicate that the extraction of Pi i resonances is insensitive to the analytic structure of the amplitude in the 
region below nN threshold. Although we did not show a result, we have also fitted to the old CMB amplitude, which 
is rather different from SAID-SES for W > 1 .5 GeV, and still found that the Roper two poles are stable. 
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